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PREFACE 


Beginning  with  Rao  and  Griliches  (1969) ,  several  studies  have 
compared  the  efficiency  of  various  estimators  of  regression  models 
with  first-order  autocorrelated  errors.  The  present  report,  however, 
is  the  first  to  compare  the  efficiency  of  all  principal  estimators 
with  trended  data  and  also  to  compare  their  performance  in  hypothesis 
testing. 

This  report  is  an  outgrowth  of  work  undertaken  to  g  Me  the 
analysis  of  school  district  expenditure  behavior,  but  it  should  be 
of  interest  to  applied  econometricians  working  in  many  different  sub¬ 
stantive  fields.  Earlier  results  are  reported  in  Park  and  Mitchell 
(1978).  The  present  report  adds  iterative  methods  to  the  list  of 
estimators  compared  before,  and  provides  more  extensive  and  illumi¬ 
nating  measures  of  relative  performance.  The  Rand  Corporation  pro¬ 
vided  computer  time  and  secretarial  services  to  make  this  new  report 
possible. 

This  study  is  also  forthcoming  in  the  Journal  of  Econometrics . 


ABSTRACT 


A  Monte  Carlo  study  is  made  of  the  small  sample  properties  of 
various  estimators  of  the  linear  regression  model  with  first-order 
autocorrelated  errors.  When  independent  variables  are  trended,  esti¬ 
mators  using  T  transformed  observations  (Prais-Winsten)  are  much  more 
efficient  than  those  using  T-l  (Cochrane-Orcutt) .  The  best  of  the 
feasible  estimators  is  iterated  Prais-Winsten  using  a  sum-of-squared- 
error  minimizing  estimate  of  the  autocorrelation  coefficient  p.  None 
of  the  feasible  estimators  performs  well  in  hypothesis  testing;  all 
seriously  underestimate  standard  errors,  making  estimated  coefficients 
appear  to  be  much  more  significant  than  they  actually  are. 
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I.  INTRODUCTION  AND  SUMMARY 


Several  estimators  are  commonly  used  to  estimate  the  linear 
regression  model  with  first-order  autocorrelated  errors.  This  Monte 
Carlo  study  extends  the  investigation  of  the  small-sample  properties 
of  such  estimators,  first  undertaken  by  Rao  and  Griliches  (1969),  in 
two  major  respects:  (1)  We  provide  a  svstematic  comparison  of  the 
estimation  efficiency  of  all  principal  estimators  with  trended  data 
and  (2)  we  compare  estimator  performance  in  hypothesis  testinp. ^ 

The  major  estimators  can  be  classified  according  to  (a)  whether 
T-l  or  T  transformed  observations  are  used,  (b)  whether  the  auto¬ 
correlation  coefficient  p  is  known  or  estimated,  and,  if  estimated, 

(c)  whether  the  estimate  of  p  is  iterated.  With  trended  data  we  find 
that  estimators  using  T-l  observations  have  very  low  efficiency,  often 
less  than  that  of  ordinary  least  squares  (OLS) ,  regardless  of  whether 
o  is  known  or  estimated.  For  unknown  p,  iterative  estimators  using 
all  T  observations  dominate  OLS  and  are  somewhat  more  efficient  than 
two-step  estimators.  We  find  that  the  iterated  Prais-Winsten  esti¬ 
mator  using  the  sum-of-squares  minimizing  p  estimate  performs  mar¬ 
ginally  better  than  the  full  maximum  likelihood  estimator. 

For  the  empirical  researcher,  reliable  hypothesis  testing  proce¬ 
dures  are  as  important  as  efficient  coefficient  estimates.  Perhaps 
the  most  serious  deficiency  of  OLS  in  the  presence  of  autocorrelation 
is  not  inefficiency  but  bias  in  its  estimated  standard  errors — a 
bias  that  in  many  situations  will  make  the  estimated  coefficients 
appear  to  be  much  more  significant  than  they  actually  are.  Unfor¬ 
tunately,  our  results  show  that  in  this  regard  the  preferred  esti¬ 
mators,  though  substantially  better  than  OLS,  can  still  be  seriously 
misleading. 


There  has  been  a  good  deal  of  recent  work  on  particular  aspects 
of  estimating  the  autocorrelated  error  model,  but  none  treats  all  of 
the  principal  estimators  and  none  deals  with  hypothesis  testing.  See 
the  discussion  below  (p.  8). 
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The  estimation  problem  and  the  estimators  that  we  consider  are 
described  in  the  next  section.  For  short  time  series  (T=20) ,  Section 
III  compares  the  efficiency  of  the  various  estimators  and  Section  IV 
describes  their  performance  in  hypothesis  testing.  Results  for 
longer  time  series  (T=50)  are  presented  in  Section  V.  Section  VI 
is  a  concise  list  of  recommendations  based  on  our  results. 
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II.  THE  ESTIMATION  PROBLEM 


THE  MODEL 

The  commonly  encountered  econometric  model  is: 


yt  =  Xt3  +  Ut  ’ 

ut  =  Put_ i  ^  =  2,  T;  (1) 

where  |p|  <  1,  E(e)  =  0,  E(ee')  =  I.  In  general,  will  include 
a  1  for  the  constant  term;  that  is. 


x 


t 


[1, 


2,  t’ 


XK,t] 


(2) 


For  this  model  the  T  x  T  covariance  matrix  of  the  error  vector  is 


E(uu' )  =  a2  V  =  a2 
u  u 


P 

1 

P 


2 

P 

P 

1 


T-l  T-2  T-3 

P  P  P 


2  2  2 
where  =  a^/Cl-p  ). 

If  p  is  known,  the  AITKEN  estimator 


b  =  (x'v^x)'1  X'V  \ 


T-l 

p 

T-2 

P 

T-3 

P 


(3) 


(4) 


is  best  linear  unbiased.  Computationally,  it  is  convenient  to 
decompose  V  ^  =  [l/(l-p ^)]R'R, 


I 
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where 


and  calculate 

b  ■=  (X'R'RXT1  X'R'Ry  (6) 

as  an  OLS  regression  of  the  transformed  variables  y*  =  Ry  on  X *  =  RX. 
THE  ESTIMATORS 

The  estimators  that  we  consider  may  all  be  thought  of  as  variants 
on  the  AITKEN  estimator.  As  shown  in  Fig.  1,  some  use  T  transformed 
observations,  and  some  use  T— 1 ;  they  also  use  different  values  for  p. 

Estimated 
Autocorrelation 
Coefficient  (p) 

Zero 
True  p 

Sutn-of- squared- 
error  minimizing 

Likelihood 
maximizing 

Fig.  1 — Estimators  considered  in  this  report. 


Number  of  Observations 
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OLS 

TRUECO 

AITKEN 

2SC0 

2SPW 

ITERCO 

ITERPW 

— 

BM 
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It  is  common  to  omit  the  first  row  in  the  transformation  matrix 
R.  ^  We  denote  the  reduced  matrix  by  S.  Then  the  transformed  vari¬ 
ables  if*  ~  Su  and  X*  =  SX  are  the  T-l  weighted  first  differences 


p  y 


t-i 


= 
'  t 


-p. 


X2,t  px2,t-l’ 


XK, t_pXK, t-l 


(7) 


This  is  the  transformation  first  proposed  by  Cochrane  and  Orcutt  (CO, 
1949)  . 

Alternatively,  Prais  and  Winsten  (PW,  1954)  recommend  retaining 
the  first  row  of  R,  in  which  case  one  has  T  transformed  observations, 
including  in  addition  to  (7)  the  transformed  first  observation 


x*  =  [  Vl-p2, 


2 

P  x 


2,1’ 


(8) 


If  the  true  value  of  p  were  known,  its  use  in  R  with  all  T  obser¬ 
vations  would  yield  the  AITKEN  estimator.  Using  true  p  and  T-l  obser¬ 
vations  would  give  what  we  call  the  TRUECO  (true  p,  Cochrane-Orcutt) 
estimator. 

In  practice,  p  is  almost  never  known.  It  is  common  to  substitute 

a  consistent  estimate  p  based  on  the  residuals  u  from  a  first-stage 

OLS  regression  using  untransformed  variables.  The  estimators  we  use 

minimize  the  sum  of  squared  errors  for  the  transformed  regression, 

2 

conditional  on  given  estimates  of  g.  For  the  CO  transformation,  the 
estimator  is 


This  is  done,  for  example,  in  the  widely  used  regression  package 
ISP;  the  TSP  procedure  CORC  is  the  same  as  what  we  call  ITERCO. 

2 

We  are  grateful  to  dames  MacKinnon  for  suggesting  that  we  use  the 
sum-of-squared-errors  minimizing  estimate  of  p.  See  below  (p.  6)  for  a 
discussion  of  alternative  estimators. 
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T  T-l  2 

>C0  =  t=2  ^ 


and  for  the  PW  transformation  it  is 

T  T-l  2 

ppw  =  tf2  Gt  at-i7tf2  Gt 


(9a) 


(9b) 


Using  p  in  S  produces  what  we  call  2SC0  (two-stage  Cochrane-Orcutt) 
estimates;  using  pp^  in  R  produces  2SPW  (two-stage  Prais-Winsten) 
estimates . 

Iterative  estimates  based  on  estimated  p  are  obtained  as  follows 
(a)  Use  the  second  stage  estimate  of  B  to  calculate  new  residuals 
u  =  y  -  Xb.  (b)  Use  these  to  calculate  a  new  estimate  of  p.  (c)  Use 
the  new  p  in  S  (or  R)  to  calculate  a  new  estimate  of  B.  (d)  Repeat 
these  steps  until  successive  estimates  of  p  differ  by  less  than  +6. 

We  set  6  =  .00001  and  call  the  resulting  estimators  ITERCO  (iterated 
Cochrane-Orcutt)  using  S  and  ITERPW  (iterated  Prais-Winsten)  using  R. 

There  is  some  chance  that  P  estimated  according  to  (9a)  or  (9b) 
will  take  on  inadmissible  values  (|  p  |  >  1) •  When  p  _>  1,  we  reset 
it  to  .99999  (=  1  -  6);  p  -  1  becomes  -.99999. 

Finally,  Beach  and  MacKinnon  (BM,  1978)  proposed  a  full  maximum 

likelihood  estimator.  Because  the  log  likelihood  function  includes  _ 

2 

the  term  .51og(l  -  p  ),  the  estimated  p  is  bounded  away  from  +  1,  so 
that  inadmissible  values  of  p  do  not  occur.  Computationally,  the  BM 
procedure  is  the  same  as  ITERPW,  except  that  a  different  estimate  of 
p  is  used;  the  BM  estimate  of  p  maximizes  the  likelihood  function 
conditional  on  estimated  8. 


ALTERNATIVE  ESTIMATORS  OF  p 

Many  previous  Monte  Carlo  studies  have  used  the  following  esti¬ 
mate  of  p  : 

T  T 

P  =  £  u  u  ,/  E  u^  .  (10) 

t=2  C  t=2 


< 

•1 
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This  estimator  is  consistent ,  but  unlike  (9a)  or  (9b)  it  does  not 
minimize  the  sum-of-squared-errors  for  either  CO  or  PW. 

In  preparing  this  report  we  found  (10)  to  be  inferior.  Using 
(9b)  rather  than  (10)  in  2SPW  and  ITERPW  reduces  the  well-known 
downward  bias  in  estimated  p,  and  results  in  slightly  smaller  root 
mean  squared  errors  for  both  p  and  b  in  almost  all  cases.  Using  (9a) 
rather  than  (10)  makes  little  difference  in  2SC0,  but  has  a  large 
effect  in  ITERCO  on  the  number  of  times  p  sticks  at  the  boundary 
value  .99999.  With  (10),  p  =  .99999  in  a  large  fraction  of  the 
experiments  in  which  true  p  _>  .8 — over  50  percent  in  one  case.  With 
(9a),  the  fraction  never  exceeded  4  percent,  and  it  was  usually  much 
smaller  than  that.^  Since  boundary  estimates  of  p  result  in  very  bad 
estimates  of  the  intercept  coefficient  6^  (for  reasons  discussed  in 
Section  III  below),  (9a)  is  decidedly  preferable  to  (10)  in  ITERCO. 

Other  consistent  estimators  of  p  have  been  pioposed.  Theil 
(1971,  p.  254)  suggests  incorporating  a  degrees-of-freedom  correction 
that  would  yield  estimates  smaller  in  absolute  value  than  (9a)  and 
(9b).  In  light  of  the  downward  bias  in  (9a)  and  (9b)  as  they  stand, 
this  correction  does  not  seem  desirable. 

Durbin  (1960)  proposes  running  an  auxiliary  regression  to  esti¬ 
mate  p.  Using  (9a)  in  ITERCO  and  (9b)  in  either  2SPW  or  ITERPW 

almost  always  resulted  in  p  estimates  that  were  less  biased  and  had 

2 

smaller  mean  squared  errors  than  the  Durbin  estimates. 

On  balance,  the  sum-of-squared-errors  minimizing  p  estimators 
(9a)  and  (9b)  appear  to  be  better  than  any  of  the  commonly  used 
alternatives. 

^See  Appendix  Table  A. 8.  The  reason  for  the  difference  is  that 
u^  tends  to  be  larger  than  u|  in  ITERCO,  because  of  the  relatively 
small  weight  given  the  first  observation  in  the  CO  transformed 
regression. 

2 

Compare  Appendix  Table  A. 6  with  Table  2  in  Park  and  Mitchell 
(1978,  p.  12). 


INDEPENDENT  VARIABLES 


We  analyze  the  relative  performance  of  the  estimators  in  Fig.  1 
using  three  independent  variables.  They  are  one  artificial  trended 
series : 

(a)  xt  =  [1,  t]; 
one  real  trended  series: 

(b)  xt  =  [1,  GNPt],  the  annual  U.S.  gross  national  product 
in  constant  dollars  beginning  in  1950; 

and,  for  comparison,  one  real  untrended  series: 

(c)  xt  =  [1,  CAP  ] ,  the  annual  U.S.  manufacturing  capacity 
utilization  rate,  also  beginning  in  1950. 

We  work  chiefly  with  20  observations  (T=20),  a  sample  size  representa¬ 
tive  of  many  time  series  studies.  In  Section  V  we  discuss  the  differ¬ 
ences  found  when  50  (quarterly)  observations  are  available. 

Strictly  speaking,  our  results  are  all  conditional  on  the  par¬ 
ticular  X  matrices  we  have  used.  But  we  believe  that  our  findings  are 
generally  applicable  for  trended  independent  variables,  because  the 
results  are  very  much  the  same  for  the  artificial  and  the  real  trended 
series.  Each  type  of  series  answers  questions  left  open  by  the  other. 
The  artificial  series  clearly  establishes  the  effect  of  pure  trends, 
but  leaves  open  the  question  of  whether  the  results  would  hold  for  the 

quirks  present  in  real-world  data — a  question  that  the  GNP  results 

2 

answer  in  the  affirmative. 

OTHER  RECENT  WORK 

Although  several  recent  papers  have  discussed  aspects  of  esti¬ 
mating  the  autocorrelated  error  model,  this  is  the  first  to  provide 
a  unified  investigation  of  all  of  the  major  estimators  with  trended 
data.  Furthermore,  with  the  exception  of  Park  and  Mitchell  (1978), 
none  of  the  previous  work  has  taken  up  the  question  of  hypothesis 
testing. 


^Maeshiro  (1976)  also  used  (a)  and  (b).  Instead  of  (c),  he  used 
quarterly  capacity  utilization  starting  with  1948. 

2 

Furthermore,  Monte  Carlo  experiments  with  quite  different  arti¬ 
ficial  time  series  yielded  similar  results;  see  Park  and  Mitchell  (1979). 
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In  a  pair  of  articles,  Maeshiro  compares  the  efficiency  of  esti¬ 
mators  using  knoun  p.  In  (1976),  he  shows  that  TRUECO  is  less  effic¬ 
ient  than  OLS  with  trended  data,  and  in  (1979)  he  demonstrates  that 
AITKEN  is  often  substantially  more  efficient  than  TRUECO.  Beach  and 
MacKinnon  (1978)  compare  their  proposed  full  maximum  likelihood  esti¬ 
mator,  BM,  with  ITERCO,  and  find  BM  to  be  more  efficient,  especially 
when  the  data  are  trended.  Harvey  and  McAvinchey,  in  an  unpublished 
paper  (1978),  make  efficiency  comparisons  of  most  of  the  major  esti¬ 
mators  applied  to  both  trended  and  untrended  data.  They  do  not, 
however,  consider  ITERPW,  the  estimator  that  we  find  to  be  the  best 
performer.  Park  and  Mitchell  (1978)  do  not  consider  any  iterative 
procedures.  Using  untrended  data,  Spitzer  (1979)  revisits  the  esti¬ 
mators  investigated  by  Rao  and  Griliches  (1969) — these  do  not  include 
ITERCO  and  ITERPW— and  adds  BM. 


III.  EFFICIENCIES  OF  ESTIMATORS 


EXACT  THEORETICAL  EFFICIENCIES 

We  can  make  two  of  our  efficiency  comparisons  using  exact 
formulas.  For  the  case  of  known  values  of  p,  the  exact  variances 
of  the  OLS ,  TRUECO,  and  AITKEN  estimators  are  given  by  the  formulas: 


var (2>  )  =  o2(X'X)~1X'VX(X'X)~ 1 

OLS  u 


(id 


Var<tTROEOO)  ' 


Vat(iAITKE»)  '  „ hx’R'RX)-1 


(12) 

(13) 


For  these  estimators  we  define  relative  efficiency  as  the  square 
root  of  the  ratio  of  the  variances  of  the  estimators  being  compared. 
For  example : 


EFF(bl,CO)  [var(bl,OLS)/var(bl,TRUECO)] 

This  definition  is  in  accord  with  comparisons  of  standard  errors  or 
t-ratios  commonly  used  by  applied  researchers;  to  use  the  ratio  it¬ 
self  would  make  the  difference  between  estimators  appear  larger  than 
they  "really"  are. 


EXPERIMENTAL  EFFICIENCIES 

We  used  Monte  Carlo  simulation  to  assess  the  relative  efficien¬ 
cies  for  the  other  five  estimators — 2SC0,  ITERCO,  2SPW,  1TERPW,  and 
BM.  For  each  of  the  independent  variables  x  and  for  each  value  of 
p  =  -.8,  .0,  .4,  .8,  .9,  .98,  we  generated  1000  samples  using  model 
(1)  with  B  =  [1,  1].  A  value  was  generated  by  drawing  a 
random  from  N(0,  1)  and  dividing  by  Vl-p.  Successive  values 
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of  et  drawn  from  N(0,  1)  were  used  to  calculate  ut  =  put  +  et> 
and  hence  y  =  a^B  +  ut*  t*len  applied  each  estimation  method  and 
averaged  the  squared  errors  of  the  estimated  coefficients  over  the 
1000  samples.^-  For  these  estimators  we  define  relative  efficiency 
as  the  ratio  of  the  root  mean  squared  errors  of  the  estimators 
being  compared.  For  example: 

EFF(bl,2SC0>  '  mSE(bl,0LS)/RMSE<bl  ,2SC0^  ’ 


where 


RMSE(b^)  =  [^^(bj-B^/lOOOT5  . 

EFFICIENCY  OF  ESTIMATORS  THAT  USE  T-l  TRANSFORMED  OBSERVATIONS 

Table  1  shows  the  efficiency,  relative  to  OLS ,  of  the  three 

estimators  that  use  T-l  observations.  Here  we  focus  on  the  results 
2 

for  positive  p.  For  trended  variables  all  three  estimators  are  less 
efficient  than  OLS  in  almost  all  of  the  cases  tabulated. 

For  xt  =  [1,  t],  TRUECO  has  extremely  low  efficiency  as  p 
approaches  1.  This  poor  performance  is  the  result  of  collinearity ; 
because  the  transformed  variable  x*  =  x„  -px_  .  approaches  the 

w  }  1  t  j  L  ^  f  L“ 1 

same  value  for  all  t,  the  T-l  vector  cc*  becomes  collinear  with  the 
transformed  vector  for  the  constant  term  xf  =  1-p ,  a  situation  we 11- 
known  for  producing  inefficient  estimates.  More  generally,  the  CO 
transformation  of  any  linearly  trended  variable  (using  p  >  0)  pro¬ 
duces  observations  that  are  more  nearly  constant  than  are  the 

^The  calculations  were  done  in  double  precision  on  an  IBM  370- 
158  using  regression  analysis  subroutines  from  the  STATLIB  statistical 
package . 

2 

Results  for  p  *  -.8  and  p  *  .0  are  included  in  the  Appendix 
tables. 

ii 


EFFICIENCY,  RELATIVE  TO  OLS,  OF  ESTIMATORS 
THAT  USE  T-l  TRANSFORMED  OBSERVATIONS 
(T=20) 


Independent 

Variable  Estimator 

P 

.4 

.8 

.9 

.98 

bl 

b2 

bl 

b2 

bl 

b2 

bl 

b2 

t 

TRUECO 

.81 

.86 

.50 

.62 

.29 

.42 

.04 

.11 

2  SCO 

.81 

.86 

.64 

.77 

.31 

.62 

.  66 

.74 

ITERCO 

.80 

.85 

.51 

.69 

.27 

.56 

.54 

.64 

GNP 

TRUECO 

.88 

.91 

.71 

.81 

.57 

.75 

.29 

.71 

2SC0 

.91 

.93 

.84 

.91 

.87 

.95 

.95 

1.03 

ITERCO 

.73 

.85 

.59 

.80 

.51 

.83 

.52 

.88 

CAP 

TRUECO 

1.10 

1.10 

1.85 

1.83 

2.10 

2.19 

1.04 

2.51 

2SC0 

1.05 

1.04 

.01 

1.41 

.00 

1.65 

.00 

1.83 

ITERCO 

1.03 

1.03 

.01 

1.65 

.00 

2.03 

.00 

2.27 

Note;  Exact  theoretical  relative  efficiency  for  TRUECO;  experi¬ 
mental  relative  efficiency  for  2SC0  and  ITERCO. 
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untransf ormed  values,  and  hence  more  nearly  collinear  with  the  con¬ 
stant  vector.^ 

When  p  is  large,  2SC0  using  trended  data  performs  better  than 

2 

TRUECO.  However,  it  is  less  efficient  than  OLS  in  almost  all  cases. 
Iterating  on  p  makes  matters  worse.  For  trended  variables,  ITERCO  is 
less  efficient  than  2SC0  in  all  of  the  cases  tabulated. 

For  the  untrended  variable  CAP  the  picture  is  mixed.  All  three 
CO  slope  estimators  are  more  efficient  than  OLS,  but  2SC0  and  ITERCO 
produce  very  bad  intercept  estimates  when  p  takes  on  the  boundary 
value  .99999.  This  happened  between  3  and  35  times  out  of  1000  trials 

3 

in  the  cases  tabulated,  causing  very  large  root  mean  squared  errors 
and  resulting  in  near-zero  relative  efficiency. 

Because  a  few  very  bad  estimates  can  dominate  the  experimental 
relative  efficiences.  Table  1  might  conceivably  hide  a  good  perform¬ 
ance  in  most  of  the  trials.  This  is  not  the  case  for  trended  vari¬ 
ables,  although  it  is  true  for  CAP.  Table  2  shows  the  number  of  times 
out  of  1000  that  each  estimator  came  closer  to  the  true  value  of  the 
coefficient  than  did  OLS.  For  trended  variables,  none  of  the  CO 
estimators  was  closer  than  OLS  as  much  as  half  of  the  time  in  any 
case  tabulated. 

_ 

Maeshiro  (1976)  explains  the  poor  performance  of  the  TRUECO 
estimator  as  a  result  of  reduction  of  variance  of  the  transformed  inde¬ 
pendent  variables.  In  one  sense  this  is  equivalent  to  our  "collinear- 
ity"  explanation,  because  a  low-variance  independent  variable  is  repre¬ 
sented  by  a  nearly  constant  vector,  which  is  nearly  collinear  with  the 
vector  for  the  constant  term.  But  in  another  sense  it  is  somewhat  mis¬ 
leading.  In  a  model  without  a  constant  term,  an  independent  variable 
with  low  (or  even  zero)  variance  causes  no  problems.  For  example,  the 
estimator  of  8  in  y  =  Bx  +  u  ,  E(u  )  =  0,  when  x  =  k  for  all  t  has 

9  2  ^  C  ^ 

variance  var(b)  =  az/Tk  . 

2 

This  refutes  Maeshiro's  (1976)  conjecture  that  "an  estimator 
utilizing  relevant  extraneous  information  (e.g.,  the  true  value  of  p) 
cannot  be  less  efficient  than  an  estimator  not  utilizing  it."  The 
reason,  paradoxically,  is  the  downward  bias  in  p  (see  Appendix  Table 
A. 7).  The  higher  the  value  of  p  used  in  the  transformation,  the  more 
collinear  are  the  transformed  variables. 

3 

See  Appendix  Table  A. 8. 

4 

The  proportion  is  significantly  less  than  one-half  on  a  binominal 
test  at  the  .05  significance  level  in  all  but  two  cases. 
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Table  2 


NUMBER  OF  TIMES  IN  1000  TRIALS  THAT  ESTIMATORS  THAT 
USE  T-l  TRANSFORMED  OBSERVATIONS  BEAT  OLS 
(T=20) 


P 


Independen  t 
Var i able 

Estimator 

.4 

.8 

.9 

.98 

bl 

b2 

bl 

b2 

bl 

b2 

bl 

b2 

t 

TRUECO 

371 

394 

236 

316 

151 

222 

24 

74 

2SCO 

382 

413 

325 

378 

324 

349 

352 

397 

ITERCO 

381 

410 

321 

372 

318 

340 

347 

383 

GNP^ 

TRUECO 

443 

461 

382 

441 

324 

398 

192 

374 

2SC0 

443 

457 

421 

468 

433 

454 

438 

495 

ITERCO 

441 

453 

405 

451 

395 

423 

395 

450 

CAP 

TRUECO 

568 

563 

726 

731 

728 

713 

523 

781 

2SC0 

547 

534 

737 

737 

742 

741 

721 

793 

ITERCO 

532 

521 

722 

725 

720 

732 

703 

786 

Note:  Counts  greater  than  531  or  smaller  than  469  are  significantly 
different  from  500  at  the  .05  level. 


-15- 


The  lesson  is  clear:  Avoid  using  any  form  of  the  Cochrane- 
Orcutt  estimator. 

EFFICIENCY  OF  ESTIMATORS  THAT  USE  T  TRANSFORMED  OBSERVATIONS 

Turning  now  to  estimators  that  use  all  T  transformed  observa¬ 
tions,  we  see  in  Table  3  that  they  all  provide  more  efficient  esti¬ 
mates  than  does  OLS.  By  retaining  the  differentially  weighted  first 
observation,  these  estimators  break  the  collinearity  that  plagued 
the  CO  estimators. 

For  trended  data,  the  AITKEN  estimator  provides  respectable, 
although  not  spectacular,  efficiency  improvements  over  OLS,  ranging 
up  to  26  percent  in  the  cases  tabulated.  The  efficiency  gain  is 
highest  in  just  those  cases  where  CO  performs  worst,  that  is,  for 
p  >_  .8.  The  other  methods,  using  estimated  rather  than  true  p,  pre¬ 
serve  about  half  of  the  AITKEN  efficiency  improvement.  For  untrended 
data,  the  efficiency  gains  are  larger,  and  more  of  the  gain  is 
retained  when  p  must  be  estimated. 

Iteration  helps.  ITERPW  is  slightly  more  efficient  than  2SPW 
using  trended  data,  and  substantially  more  efficient  with  untrended 
data.  The  BM  estimator  has  virtually  the  same  efficiency  as  ITERPW. 

Using  ITERPW  as  a  standard  of  comparison,  we  show  in  Table  4  the 
number  of  times  in  1000  trials  that  the  other  estimators  using  T 
transformed  observations  come  closer  to  the  true  coefficient  value. 

In  almost  all  cases,  ITERPW  outperformed  the  other  estimators  that 
use  estimated  p. 


i 
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Table  3 

EFFICIENCY,  RELATIVE  TO  OLS,  OF  ESTIMATORS  THAT 
USE  T  TRANSFORMED  OBSERVATIONS 
(T=20) 


Independent 
Var i able 

Estimator 

_ e _ 1 

.4 

.8 

.9 

.98 

bl 

b2 

bl 

b2 

bl 

b2 

bl 

b2 

t 

AITKEN 

1.02 

J  .02 

l.ub 

1.09 

1.08 

1.10 

1.03 

1.08 

2SPW 

1.01 

1.01 

1.05 

1.06 

1.05 

1.08 

1.02 

1.05 

ITERPW 

1.01 

1.01 

1.06 

1.07 

1.05 

1.08 

1.02 

1.05 

BM 

1.01 

1.01 

1.05 

1.06 

1.05 

1.07 

1.02 

1.04 

GNPt 

AITKEN 

1.02 

1.02 

1.13 

1.14 

1.16 

1.20 

1.13 

1.26 

2SPW 

1.01 

1.01 

1.08 

1.09 

1.10 

1.13 

1.06 

1.12 

ITERPW 

1.00 

1.01 

1.08 

1.09 

1.11 

1.14 

1.07 

1.14 

BM 

1.01 

1.01 

1.08 

1.09 

1.10 

1.13 

1.06 

1.12 

CAP 

AITKEN 

1.14 

1.14 

1.85 

1.86 

2.15 

2.21 

1.95 

2.52 

L 

2SPW 

1.06 

1.06 

1.38 

1.38 

1.57 

1.60 

1.54 

1.75 

ITERPW 

1.04 

1.05 

1.61 

1.61 

1.94 

2.00 

1.77 

2.15 

BM 

1.05 

1.06 

1.58 

1.58 

1.92 

— L-92 

1.75 

2.12 

Note :  Exact  theoretical  relative  efficiency  for  AITKEN;  experimental 

relative  efficiency  for  2SPW,  I1ERPW,  and  BM. 


I 
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Table  4 


NUMBER  OF  TIMES  IN  1000  TRIALS  THAT  ESTIMATORS  THAT 
USE  T  TRANSFORMED  OBSERVATIONS  BEAT  ITERPW 
(T=20) 


p 


Independent 

Variable  Estimator 

.4 

.8 

.9 

.98 

bl 

b2 

bl 

b2 

bl 

b2 

bl 

b2 

t 

AITKEN 

517 

503 

534 

537 

517 

546 

563 

542 

2SPW 

493 

514 

449 

445 

430 

420 

458 

439 

BM 

489 

510 

446 

445 

450 

406 

433 

432 

GNP 

AITKEN 

532 

532 

527 

527 

564 

591 

568 

550 

2SPW 

502 

499 

433 

428 

391 

383 

419 

377 

BM 

499 

494 

438 

433 

382 

380 

409 

364 

CAP 

AITKEN 

561 

565 

561 

564 

523 

531 

539 

538 

2SPW 

500 

500 

398 

390 

399 

423 

434 

417 

BM 

512 

514 

430 

419 

436 

449 

427 

423 

Note:  Counts  greater  than  531  or  smaller  than  469  are  significantly 
different  from  500  at  the  .05  level. 
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IV.  HYPOTHESIS  TESTING 


In  our  opinion,  the  most  troublesome  characteristic  of  OLS  when 

V  I  is  not  the  loss  in  efficiency,  but  the  bias  in  the  statistic 
-2  -1 

a^(X'X)  ,  which  is  conventionally  used  as  an  estimator  of  the  co- 

variance  matrix  of  the  estimated  coefficients.  How  serious  this 
bias  can  be  is  illustrated  by  the  experimental  results  shown  in 
Table  5.  We  focus  on  the  methods  that  use  all  T  transformed  obser¬ 
vations,  since  they  are  always  more  efficient  than  methods  using  T-l 
observations.  OLS  seriously  underestimates  standard  errors  for  all 
cases  shown.  For  example,  for  trended  data  and  p  >^  . 8 ,  when  one 
applies  a  two-tailed  test  at  the  .05  level,  the  underestimate  is 
large  enough  to  lead  one  to  judge  an  estimated  coefficient  to  be 
significantly  different  from  its  true  va’ue  45  to  85  percent  of  the 
time. 

The  AITKEN  estimator  provides  unbiased  variance  estimates,  of 
course.  Unfortunately,  the  procedures  using  estimated  p  do  not  do 
as  well,  although  they  do  improve  on  OLS.^  ITERPW  is  the  best  of  the 
lot.  Still,  for  p  ^  .8  and  trended  data,  it  would  result  in  rejec¬ 
tion  of  a  correct  null  hypothesis  at  least  25  percent  of  the  time. 

For  the  untrended  variable  CAP,  the  results  are  qualitatively 
similar,  but  the  biases  are  less  severe  than  in  the  case  of  trended 
variables. 

On  hypothesis  testing  grounds,  ITERPW  appears  to  be  superior  to 
both  2SPW  and  BM.2 

^The  covariance  matrix  of  the  estimated  coefficients  is  esti- 
mated  directly  in  the  transformed  regressions  as  o£(X*'X*)  , 

where  o2  =  e'e/CT-K). 

2  ~2  »  „ 

Had  we  used  the  maximum  likelihood  estimate  o  =  e'e/T  in  the 

BM  procedure,  the  margin  of  superiority  of  ITERPW  over  BM  would  have 

been  slightly  larger. 
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TabJe  5 


NUMBER  OK  TYPE  1  ERRORS  IN  1000  TRIALS 
AT  .05  SIGNIFICANCE  LEVEL 
(T=20) 


Kst imat  or 

p 

I ndependen  t 
Variable 

.4 

.8 

.9 

.98 

bl 

b2 

bl 

b2 

h 

b2 

bl 

b2 

t 

OLS 

193 

197 

502 

490 

645 

571 

848 

709 

2SPW 

125 

132 

302 

293 

411 

340 

690 

473 

ITERPW 

124 

131 

293 

285 

401 

336 

700 

474 

BM 

126 

133 

312 

305 

433 

360 

731 

503 

CNP 

OLS 

186 

185 

457 

449 

601 

596 

730 

666 

t 

2SPW 

138 

138 

258 

251 

354 

343 

509 

413 

ITERPW 

136 

136 

254 

246 

343 

322 

486 

397 

BM 

139 

138 

261 

258 

375 

352 

534 

432 

CAP 

OLS 

147 

143 

304 

294 

341 

323 

407 

322 

t 

2SPW 

110 

106 

153 

149 

154 

144 

211 

137 

ITERPW 

115 

113 

102 

101 

90 

86 

144 

86 

BM 

112 

109 

107 

107 

98 

92 

176 

86 

Note:  Counts  greater  than  63  or  smaller  than  37  are  significantly 
different  from  50  at  the  .05  level. 
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V.  LONGER  TIME  SERIES 


In  this  section  we  investigate  how  the  results  change  when  50, 
rather  than  20,  observations  are  available  for  one  of  our  independent 
variables,  x  =  [1,  GNPt].  In  order  to  increase  the  length  of  the 
GNP  time  series  to  get  50  observations,  we  must  shift  from  annual  to 
quarterly  observations.  Because  the  quarterly  series  exhibits  short¬ 
term  fluctuations  that  are  averaged  out  in  the  annual  series,  our 
T=50  series  is  less  trended  than  the  annual  GNP  series  used  above. 

It  is,  however,  typical  of  longer  economic  time  series. 

A  larger  sample  markedly  improves  the  estimators  of  p.  For 

example,  when  true  p  =  .9,  the  mean  value  of  p........  increased  from 

1 1  LKr  W 

.59  for  T=20  to  .80  for  T=50.  The  bias,  although  still  clearly 
apparent,  is  greatly  reduced. 

Tables  6,  7,  and  8  repeat  the  information  in  Tables  1  through  5 
for  T=50.  For  the  most  part,  the  conclusions  for  T=20  also  apply  for 
the  larger  sample  size.  Estimators  using  T-l  transformed  observa¬ 
tions  are  usually  less  efficient  than  0LS .  Those  using  T  observa¬ 
tions  always  improve  on  OLS ,  and  the  margin  is  wider  for  the  larger 
sample  size.  Also,  methods  using  estimated  p  retain  more  of  the 
AITKEN  estimator's  margin  of  improvement  (reflecting  the  improved 
estimate  of  p) .  ITERPW  appears  to  be  slightly  better  than  either 
2SPW  or  BM. 

Increased  sample  size  does  nothing  to  reduce  the  bias  in  the 
OLS  estimated  standard  errors,  but  it  does  improve  the  ITERPW  esti¬ 
mates.  Nevertheless,  ITERPW  would  still  lead  to  rejection  of  a 
correct  null  hypothesis  up  to  30  percent  of  the  time  in  the  cases 
tabulated . 

^See  Appendix  Table  A. 6 . 
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Table  6 


EFFICIENCY  COMPARISONS  FOR  ESTIMATORS  THAT 
USE  T-l  TRANSFORMED  OBSERVATIONS 
(T=50) 


Independent 

Variable 

Estimator 

D 

.4 

.8 

.9 

.98 

bl 

b2 

bl  b2 

bl 

b2 

bl 

b2 

Efficiency  relative 

to  OLS 

a 

GNP 

TRUECO 

m 

.84  .87 

.87 

.94 

.92 

1.39 

2SC0 

m m 

.85  .88 

.91 

.95 

.02 

1.13 

ITERCO 

.91 

.70  .77 

.72 

.84 

.02 

1.12 

Number  of  times  in  1000 

trials 

that 

other  estimators  beat  0LSb 

GNP 

TRUECO 

m 

428 

406  430 

433 

461 

445 

599 

2SC0 

mm 

432 

401  416 

449 

477 

509 

553 

ITERCO 

WM 

430 

387  400 

399 

429 

423 

512 

£ 

Exact  theoretical  relative  efficiency  for  TRUECO;  experimental 
relative  efficiency  for  2SC0  and  ITERCO. 

^Counts  greater  than  531  or  smaller  than  469  are  significantly 
different  from  500  at  the  .05  level. 
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Table  7 


EFFICIENCY  COMPARISONS  FOR  ESTIMATORS  THAT 
USE  T  TRANSFORMED  OBSERVATIONS 
(T=50) 


Independent 
Var i able 

Estimator 

P 

, 

4 

8 

.9 

.98 

bl 

b2 

bl 

b2 

bl 

bl 

b2 

Efficiency  relative 

to  0LS‘ 

l 

GNP 

AITKEN 

1.02 

1.02 

1.19 

1.19 

1.40 

1.41 

1.78 

1.89 

2SPW 

1.02 

1.02 

1.15 

1.14 

1.28 

1.28 

1.41 

1.44 

ITERPW 

1.02 

1.02 

1.12 

1.11 

1.26 

1.26 

1.52 

1.57 

BM 

1.02 

1.02 

1.14 

1.13 

1.28 

1.29 

1.48 

1.52 

Number  of  time  in 

1000  trials  that 

other  estimators  beat 

ITERPWb 

GNP 

AITKEN 

494 

Bl 

556 

553 

561 

586 

597 

2SPW 

500 

491 

ESS 

432 

426 

415 

383 

354 

BM 

493 

484 

441 

445 

427 

415 

384 

353 

Exact  theoretical  relative  efficiency  for  AITKEN;  experimental 
relative  efficiency  for  2SPW,  ITERPW,  and  BM. 

^Counts  greater  than  531  or  smaller  than  469  are  significantly 
different  from  500  at  the  .05  level. 
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Table  8 


NUMBER  OF  TYPE  1  ERRORS  IN  1000  TRIALS  AT 
.05  SIGNIFICANCE  LEVEL 
(T=50) 


Independent 

Variable 

Estimator 

P  ' 

•A 

.8 

•9 

\o 

00 

bl 

b2 

bl 

b2 

bl 

b2 

bl 

b2 

GNP 

OLS 

209 

209 

501 

505 

636 

631 

754 

781 

2SPW 

90 

92 

143 

143 

198 

202 

377 

366 

ITERPW 

87 

92 

138 

137 

184 

191 

307 

296 

BM 

90 

92 

151 

151 

200 

208 

357 

338 

Note:  Counts  greater  than  63  or  smaller  than  37  are  significantly 
different  from  50  at  the  .05  level. 
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VI .  RECOMMENDATIONS 


Our  results  lead  us  to  offer  the  following  guidelines  to  prac¬ 
ticing  econometricians  working  with  trended  data  in  the  presence 
of  autocorrelation: 

1.  Avoid  the  Cochrane-Orcutt  estimator  (using  T-l  transformed 
observations) ;  it  is  more  complicated  than  OLS  and  often  less  ef¬ 
ficient  . 

2.  Use  the  it  native  version  of  the  Prais-Winsten  estimator 
(using  T  transformed  observations).  It  offers  efficiency  gains  over 
OLS  that  range  from  modest  to  substantial.  It  is  slightly  but 
clearly  superior  to  two-stage  Prais-Winsten.  For  trended  data  and 

a  large  autocorrelation  coefficient,  it  also  appears  to  have  a 
slight  edge,  in  small  samples,  over  the  full  maximum  likelihood 
method  proposed  by  Beach  and  MacKinnon. 

3.  Distrust  the  conventional  t-statistics.  The  OLS  standard 
errors  are  vastly  underestimated.  The  iterative  Prais-Winsten 
standard  errors  are  a  substantial  improvement,  but  still  highly 
misleading.  Because  estimated  coefficients  seem  much  more  signifi¬ 
cant  than  they  really  are,  apply  a  more  stringent  confidence  level 
for  hypothesis  testing. 
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APPENDIX  TABLES 


The  tables  in  the  text,  above,  show  results  for  selected 
estimators  for  positive  values  of  p  only.  The  following  Appendix 
tables  include  results  for  all  relevant  estimators  for  p  -  -.8,  .0, 

.4,  .8,  .9,  and  .98. 

Table  _ Title _ 

A. 1  Exact  theoretical  efficiency,  relative  to  OLS,  of  estimators 

that  use  true  p 

A. 2  Experimental  efficiency  of  various  estimators  relative  to 

OLS 

A. 3  Number  of  times  in  1000  trials  that  various  estimators  beat 

OLS 

A. 4  Number  of  times  in  1000  trials  that  various  estimators  beat 

ITERPW 

A. 5  Number  of  type  1  errors  in  1000  trials  at  .05  significance 

level 

A. 6  Performance  of  estimators  of  p 

A. 7  Number  of  times  in  1000  trials  that  various  estimators  of 

p  beat  the  ITERPW  estimator  of  p 

A. 8  Number  of  times  in  1000  trials  that  estimated  p  equals 

boundary  value 

A. 9  Number  of  iterations  and  failures  to  converge 
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Table  A. 2 

EXPERIMENTAL  EFFICIENCY  OF  VARIOUS  ESTIMATORS  RELATIVE  TO  OLS 


Independent 

Variable 

Estimator 

T=2U 

t 

TRUECO 

2SCO 

ITERCO 

AITKEN 

2SPW 

ITERPW 

BM 

GNPt 

TRUECO 

2SC0 

ITERCO 

AITKEN 

2SPW 

ITERPW 

BM 

CAPt 

TRUECO 

2  SCO 

ITERCO 

AITKEN 

2SPW 

ITERPW 

BM 

T-50 

GNPt 

TRUECO 

2  SCO 

ITERCO 

AITKEN 

2SPW 

ITERPW 

BM 

1.17 

1.21 

.90 

.92 

.79 

.84 

.51 

.63 

.28 

.41 

.04 

.11 

1.16 

1.19 

.91 

.92 

.81 

.86 

.64 

.77 

.31 

.62 

.66 

.74 

1.16 

1.19 

.91 

.92 

.80 

.85 

.51 

.69 

.27 

.56 

.54 

.64 

1.18 

1.21 

1.00 

1.00 

1.02 

1.02 

1.09 

1.09 

1.07 

1.11 

1.05 

1.07 

1.18 

1.21 

.99 

.99 

1.01 

1.01 

1.05 

1.06 

1.05 

1.08 

1.02 

1.05 

1.18 

1.21 

.99 

.99 

1.01 

1.01 

1.06 

1.07 

1.05 

1.08 

1.02 

1.05 

1.18 

1.21 

.99 

.99 

1.01 

1.01 

1.05 

1.06 

1.05 

1.07 

1.02 

1.04 

1.21 

1.21 

.94 

.95 

.90 

.94 

.72 

.82 

.58 

.75 

.28 

.71 

1.20 

1.20 

.94 

.95 

.91 

.93 

.84 

.91 

.87 

.95 

.95 

1.03 

1.20 

1.20 

.94 

.95 

.73 

.85 

.59 

.80 

.51 

.83 

.52 

.88 

1.22 

1.21 

1.00 

1.00 

1.02 

1.03 

1.12 

1.14 

1.17 

1.21 

1.15 

1.28 

1.21 

1.21 

.99 

.99 

1.01 

1.01 

1.08 

1.09 

1.10 

1.13 

1.06 

1.12 

1.21 

1.21 

.99 

.99 

1.00 

1.01 

1.08 

1.09 

1.11 

1.14 

1.07 

1.14 

1.21 

1.21 

.99 

.99 

1.01 

1.01 

1.08 

1.09 

1.10 

1.13 

1.06 

1.12 

1.70 

1.69 

.98 

.98 

1.12 

1.11 

1.91 

1.90 

2.09 

2.18 

1.03 

2.49 

1.69 

1.68 

.97 

.97 

1.05 

1.04 

.01 

1.41 

.00 

1.65 

.00 

1.83 

1.69 

1.68 

.95 

.95 

1.03 

1.03 

.01 

1.65 

.00 

2.03 

.00 

2.27 

1.70 

1.69 

1.00 

1.00 

1.15 

1.15 

1.90 

1.90 

2.13 

2.19 

1.92 

2.49 

1.70 

1.69 

.99 

.99 

1.06 

1.06 

1.38 

1.38 

1.57 

1.60 

1.54 

1.75 

1.70 

1.69 

.98 

.98 

1.04 

1.05 

1.61 

1.61 

1.94 

2.00 

1.77 

2.15 

1.70 

1.69 

.98 

.98 

1.05 

1.06 

1.58 

1.58 

1.92 

1.97 

1.75 

2.12 

1.13 

1.13 

.96 

.96 

.91 

.92 

.83 

.87 

.83 

.90 

.90 

1.33 

1.13 

1.13 

.96 

.96 

.91 

.92 

.85 

.88 

.91 

.95 

.02 

1.13 

1.13 

1.13 

.96 

.96 

.91 

.92 

.70 

.77 

.72 

.84 

.02 

1.12 

1.14 

1.15 

1.00 

1.00 

1.02 

1.02 

1.21 

1.20 

1.38 

1.39 

1.75 

1.84 

1.14 

1.15 

1.00 

1.00 

1.02 

1.02 

1.15 

1.14 

1.28 

1.28 

1.41 

1.44 

1.14 

1.15 

1.00 

1.00 

1.02 

1.02 

1.12 

1.11 

1.26 

1.26 

1.52 

1.57 

1.14 

1.15 

1.00 

1.00 

1.02 

1.02 

1.14 

1.13 

1.28 

1.29 

1.48 

1.52 

Table  A. 3 
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Table  A.  4 

NUMBER  OF  TIMES  IN  1000  TRIALS  THAT  VARIOUS  ESTIMATORS  BEAT  ITERPW 


Xod«p«n4ent 

Variable 

Eat  Isa  cor 

_ 

a 

.0 

•* 

■  a 

_ 

98 

bi 

»2 

bl 

b2 

bi 

b2 

bi 

b2 

bi 

b2 

bi 

*2 

T*20 

t 

OLS 

421 

399 

492 

492 

463 

483 

428 

411 

417 

369 

423 

409 

TRUECO 

470 

482 

412 

430 

387 

400 

242 

299 

149 

218 

24 

75 

2  SCO 

471 

468 

422 

431 

402 

414 

343 

363 

325 

319 

357 

38) 

ITERC0 

476 

499 

420 

430 

402 

413 

333 

357 

316 

313 

349 

367 

AITKEN 

503 

506 

492 

492 

517 

503 

534 

537 

517 

546 

563 

542 

2SPW 

521 

509 

509 

512 

493 

514 

449 

445 

430 

420 

458 

439 

BM 

515 

500 

504 

505 

489 

510 

446 

445 

450 

406 

433 

432 

GNP_ 

OLS 

393 

382 

516 

533 

475 

465 

394 

381 

356 

352 

381 

346 

TRUECO 

492 

481 

465 

465 

448 

454 

349 

401 

302 

357 

181 

342 

2  SCO 

500 

497 

457 

454 

435 

449 

375 

407 

364 

379 

416 

411 

ITERCO 

499 

495 

455 

453 

433 

444 

362 

392 

347 

370 

376 

382 

AITKEN 

506 

496 

516 

533 

532 

532 

527 

527 

564 

591 

568 

550 

2SPW 

520 

531 

538 

552 

502 

499 

433 

428 

394 

383 

419 

377 

BM 

511 

525 

529 

545 

499 

494 

438 

433 

382 

380 

409 

364 

CAP 

OLS 

302 

302 

508 

504 

441 

439 

278 

263 

250 

252 

274 

216 

TRUECO 

526 

523 

510 

509 

519 

514 

552 

547 

517 

519 

310 

521 

2SCO 

527 

520 

485 

491 

473 

464 

456 

464 

412 

427 

401 

436 

ITERCO 

525 

512 

469 

474 

458 

450 

528 

515 

456 

498 

435 

514 

AITKEN 

500 

496 

508 

504 

561 

565 

561 

564 

523 

531 

539 

538 

2SPW 

503 

501 

527 

531 

500 

500 

398 

390 

399 

423 

434 

417 

BM 

490 

486 

526 

531 

512 

514 

430 

419 

436 

449 

427 

423 

T-50 

GNP 

OLS 

420 

430 

514 

516 

465 

455 

380 

372 

326 

328 

292 

272 

TRUECO 

472 

474 

449 

451 

442 

442 

Bin 

354 

365 

344 

435 

2  SCO 

466 

469 

454 

451 

427 

431 

mn 

343 

354 

319 

345 

ITERCO 

467 

470 

454 

451 

425 

430 

347 

355 

315 

319 

299 

329 

AITKEN 

474 

480 

514 

516 

490 

494 

564 

556 

553 

561 

586 

597 

2SPW 

519 

526 

549 

550 

500 

491 

440 

432 

426 

415 

383 

354 

BM 

516 

522 

517 

519 

493 

484 

441 

445 

427 

415 

384 

353 

Note:  Counts  greater  than  531  or  smaller  than  469  are  significantly  different  from  500 
at  the  .05  level. 
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Table  A. 5 

NUMBER  OF  TYPE  1  ERRORS  IN  1000  TRIALS  AT  .05  SIGNIFICANCE  LEVEL 


Independent 

Variable  E« tine tor 

.0 

.* 

_ ,8 _ 

_ j _ 

_ 1 

bi 

b2 

bl 

b2 

>1 

b2 

>1 

b2 

bi 

b2 

bi 

b2 

T-20 

■ 

N 

■R 

IM 

t 

OLS 

0 

0 

i 

197 

490 

KSf 

8-*8 

709 

TRUECO 

44 

52 

El 

■Ea 

57 

54 

B 9 

mm 

56 

54 

2SCO 

42 

57 

Ba 

839 

139 

arm 

430 

397 

682 

509 

ITERCO 

43 

61 

81 

73 

131 

138 

SS3. 

421 

390 

670 

503 

AITKEN 

41 

51 

42 

44 

47 

65 

53 

41 

43 

46 

53 

64 

2SPW 

43 

62 

78 

72 

125 

132 

302 

293 

411 

340 

690 

473 

ITERPW 

49 

65 

78 

72 

124 

131 

293 

285 

401 

336 

700 

474 

BM 

41 

55 

78 

71 

126 

133 

312 

305 

433 

360 

731 

503 

GNP 

OLS 

0 

0 

63 

59 

186 

185 

457 

449 

601 

596 

730 

666 

TRUECO 

50 

53 

49 

56 

47 

40 

45 

51 

55 

47 

57 

49 

2  SCO 

47 

56 

86 

86 

131 

128 

278 

269 

394 

371 

510 

424 

ITERCO 

48 

57 

86 

87 

127 

129 

281 

271 

392 

378 

519 

426 

AITKEN 

49 

56 

63 

59 

47 

41 

51 

53 

46 

35 

55 

52 

2SPW 

51 

60 

88 

89 

138 

138 

258 

251 

354 

343 

509 

413 

ITERPW 

52 

64 

91 

89 

136 

136 

254 

246 

343 

322 

486 

397 

BM 

48 

56 

88 

88 

139 

138 

261 

258 

375 

352 

534 

432 

CAP 

OLS 

3 

3 

52 

51 

147 

143 

304 

294 

341 

323 

407 

322 

TRUECO 

53 

55 

55 

52 

43 

45 

50 

49 

42 

43 

50 

51 

2SC0 

53 

52 

75 

70 

101 

96 

147 

143 

156 

137 

241 

129 

ITERCO 

54 

55 

78 

75 

107 

100 

93 

96 

94 

78 

198 

85 

AITKEN 

47 

50 

52 

51 

43 

46 

51 

48 

41 

49 

53 

51 

2SPW 

56 

58 

74 

74 

no 

106 

153 

149 

154 

144 

211 

137 

ITERPW 

60 

59 

80 

78 

115 

113 

102 

101 

90 

86 

144 

86 

BM 

54 

56 

78 

74 

112 

109 

107 

107 

98 

92 

176 

86 

T-50 

GNP 

OLS 

0 

0 

49 

46 

209 

209 

501 

505 

636 

631 

754 

781 

TRUECO 

50 

50 

48 

45 

55 

54 

49 

44 

57 

56 

50 

59 

2SC0 

48 

52 

67 

72 

94 

99 

163 

164 

241 

246 

397 

375 

ITERCO 

48 

52 

66 

72 

94 

99 

176 

177 

253 

254 

376 

350 

AITKEN 

52 

48 

49 

46 

51 

45 

40 

45 

52 

55 

54 

59 

2SPW 

51 

52 

69 

74 

90 

92 

143 

143 

198 

202 

377 

366 

ITERPW 

52 

53 

69 

74 

87 

92 

138 

137 

184 

191 

307 

296 

BM 

50 

51 

69 

73 

90 

92 

151 

151 

200 

208 

357 

338 

Note:  Counts  greater  than  63  or  smaller  than  37  are  significantly  different  from  50  at  the  .05 
level. 


PERFORMANCE  OF  ESTIMATORS  OF 


.78  .09  -.05  .15  .33  .15  .69  .16  .78  .16  .83  .18 


Table  A. 7 


NUMBER  OF  TIMES  IN  1000  TRIALS  THAT  VARIOUS  ESTIMATORS  OF  p 
BEAT  THE  ITERPW  ESTIMATOR  OF  p 


Independent 

P 

Variable 

Estimator 

-.8 

.0 

.4 

.8 

.9 

.98 

T=20 

t 

2  SCO 

545 

1000 

375 

83 

35 

34 

ITERCO 

557 

821 

470 

236 

201 

226 

2SPW 

588 

1000 

389 

83 

76 

147 

BM 

557 

1000 

381 

92 

39 

14 

GNP 

2SCO 

517 

1000 

372 

76 

35 

18 

ITERCO 

530 

808 

432 

246 

230 

271 

2SPW 

566 

1000 

381 

96 

101 

114 

BM 

531 

1000 

382 

97 

52 

13 

CAP 

2SCO 

515 

1000 

427 

174 

108 

77 

ITERCO 

564 

728 

545 

372 

268 

236 

2SPW 

568 

999 

437 

200 

131 

174 

BM 

566 

1000 

454 

276 

195 

64 

T=50 

GNP 

2SCO 

511 

1000 

342 

180 

107 

37 

L 

ITERCO 

516 

760 

382 

318 

287 

298 

2SPW 

541 

1000 

351 

219 

194 

198 

BM 

512 

1000 

349 

197 

133 

39 

Note:  Counts  greater  than 

531  or 

smaller  than 

469 

are 

significantly  different 

from 

500  at 

the  .05 

level 

. 

Table  A. 8 

NUMBER  OF  TIMES  IN  1000  TRIALS  THAT  ESTIMATED  p 
EQUALS  BOUNDARY  VALUE 


Independent 

Variable 

_  _  _ e _ _ 

Estimator 

-.8 

.4 

o 

.9 

.98 

T-20 

t 

2SC0 

3 

0 

0 

0 

0 

0 

ITERCO 

11 

0 

0 

0 

0 

0 

2SPW 

14 

0 

0 

2 

5 

13 

ITERPW 

58 

0 

0 

0 

0 

0 

BM 

0 

0 

0 

0 

0 

0 

GNP 

2SCO 

0 

0 

0 

0 

0 

0 

ITERCO 

7 

0 

0 

0 

0 

0 

2SPW 

12 

0 

0 

3 

8 

14 

ITERPW 

53 

0 

0 

0 

1 

6 

BM 

0 

0 

0 

0 

0 

0 

CAP 

2SC0 

4 

0 

0 

3 

10 

35 

t 

ITERCO 

8 

0 

0 

3 

10 

35 

2SPW 

27 

0 

0 

8 

29 

95 

ITERPW 

51 

0 

0 

10 

43 

136 

BM 

0 

0 

0 

0 

0 

0 

T=50 

GNP 

2  SCO 

0 

0 

0 

0 

0 

3 

ITERCO 

0 

0 

0 

0 

0 

0 

2SPW 

1 

0 

0 

0 

2 

16 

ITERPW 

1 

0 

0 

0 

7 

37 

BM 

0 

0 

0 

0 

0 

0 

34 


o 

o 

O 

o 

CO 

o 

O 

vO 

*H 

o 

CO 

CN 

rH 

z 

nO 

00 

vO 

nO 

<r 

r- 

<0 

CN 

00 

O 

00 

rH 

CN 

00 

Q 

• 

• 

• 

• 

• 

• 

• 

• 

a 

a 

a 

on 

to 

CN 

rH 

»H 

CO 

CO 

CN 

CN 

<r 

CO 

CN 

m 

•0 

* 

rH 

i—t 

rH 

rH 

p- 

m 

CN 

00 

H- 

ON 

o 

CN 

<0 

oo 

^0 

00 

a 

O 

m 

in 

rH 

vO 

m 

CO 

NO 

NO 

rH 

p>» 

NO 

rH 

rH 

iH 

rH 

u 

00 

o 

o 

NT 

o 

o 

ON 

O 

iH 

NO 

rH 

rH 

!« 

co 

m 

CO 

rH 

n£> 

no 

CO 

ON 

0N. 

ON 

CN 

*0- 

vO 

00 

ON 

Q 

• 

• 

• 

• 

• 

• 

a 

a 

a 

• 

• 

CO 

O 

rH 

rH 

rH 

CN 

rH 

o 

CO 

CO 

00 

CO 

CO 

rH 

rH 

rH 

•<0 

CO 

o 

o 

*H 

m 

lO 

p- 

ON 

ON 

-0" 

• 

• 

• 

• 

• 

■ 

• 

a 

• 

a 

a 

a 

ON 

m 

m 

o 

nO 

in 

o 

NO 

NO 

r>. 

in 

m 

rH 

rH 

o 

o 

o 

CO 

O 

o 

n- 

CN 

rH 

p^ 

o 

O 

z 

iH 

CN 

00 

Q 

<r 

vO 

<r 

00 

CN 

<0- 

CN 

CN 

co 

p^ 

CO 

• 

• 

• 

• 

• 

• 

• 

a 

a 

a 

• 

a 

p^ 

iH 

rH 

00 

CN 

rH 

NO 

<r 

NO 

CN 

rH 

CN 

rH 

ON 

r*. 

<r 

o 

m 

o 

On 

o 

00 

NO 

S 

• 

• 

• 

• 

■ 

• 

a 

a 

a 

a 

a 

• 

r>. 

m 

<r 

m 

m 

oo 

1^. 

NO 

NO 

<0 

<0 

u 

o 

o 

o 

<r 

o 

o 

CN 

CN 

O 

o 

O 

o 

z 

N0 

Q 

ON 

CN 

o 

m 

CN 

NO 

rH 

m 

r*** 

P^ 

to 

rH 

rH 

rH 

co 

rH 

rH 

>0- 

>0- 

CO 

nO 

CO 

CN 

nO 

*0- 

CN 

>0- 

CN 

o 

in 

m 

in 

a 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

"0- 

<0 

<r 

-0- 

NO 

NO 

NO 

co 

m 

CO 

U 

O 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

55 

o 

Q 

ON 

oo 

00 

00 

ON 

oo 

ON 

<0 

CO 

*4- 

>0* 

<r 

CO 

' 

’ 

CN 

CN 

CN 

ON 

00 

00 

On 

ON 

00 

O 

00 

P-* 

CN 

«H 

a 

* 

• 

• 

• 

• 

a 

• 

a 

• 

a 

CO 

CO 

CO 

co 

CO 

CO 

m 

^0- 

>0- 

CO 

CO 

CO 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

00 

Q 

vO 

n£> 

NO 

r*. 

NO 

NO 

NO 

m 

<0 

NT 

' 

CO 

00 

00 

00 

00 

ON 

00 

rH 

CN 

CN 

a 

• 

■ 

• 

• 

• 

• 

» 

a 

a 

• 

CO 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

u 

o 

•u 

2 

o 

s 

o 

re 

o 

3 

o 

3 

B 

u 

Oj 

u 

S3 

u 

o 

o-i 

•H 

Pi 

3 

oi 

Pi 

3 

Pi 

3 

U 

W 

w 

w 

W 

w 

W 

w 

w 

05 

H 

H 

a 

H 

H 

a 

H 

H 

a 

H 

H 

a 

W 

M 

M 

CQ 

M 

H 

PQ 

M 

M 

PQ 

M 

H 

CQ 

4J 

a 

05 

05 

*T3 

*H 

c 

aO 

05 

tO 

■ 

J 

. 

J 

■ 

■J 

a  -h 

col 

PH 

Si 

,(U 

05 

z 

oiz 

O 

05 

C 

o 

OO 

•H 

c 

H 

•H 

CO 

aC 

M 

o 

05 

cO 

U 

05 

■H 

u 

IH 

05 

O 

u 

o 

H 

IH 

05 

05 

aO 

aO 

B 

D 

05 

C 

00 

u 

05 

05 

aC 

> 

u 

c 

o 

4-4 

o 

O 

o 

c 

u 

o 

*H 

XJ 

05 

«0 

rH 

•H 

♦H 

> 

CO 

0) 

IH 

TJ 

CO 

T3 

CO 

U 

05 

CO 

o 

TO 

o 

a 

cO 

a 

4J 

U5 

c 

o 

CL5 

•H 

x. 

4J 

U 

CO 

u 

CO 

05 

fH 
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